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ON THE SEASONAL AND DAY-TO-DAY VARIATIONS 
IN THE F2 REGION 


H. G. Mayr* and K. K. Mahajan* 

Aeronomy Branch, NASA Goddard Space Flight Center 
Greenbelt, Maryland 


ABSTRACT 

Thomson scatter measurements at Arecibo (18.4°N, 67 °W) indicate appre- 
ciable seasonal and day-to-day changes in the noontime F-region electron con- 
centration, with a strong positive correlation between the height of the F-layer 
peak and the peak electron concentration. This correlation is theoretically re- 
produced by changing, first the neutral particle concentration and later the 
neutral wind velocity. The changes required in the neutral concentration are 
much too high to be explained on the basis of the available information on the 
atmospheric density from the satellite drag data. Horizontal winds directed 
polewards during winter and equatorialwards during summer, however, can 
reasonably explain the observed variations and are consistent in terms of the 
expected seasonal change in the position of the pressure (density) bulge with 
respect to Arecibo. 

*NRC-NASA postdoctoral research associate. 
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ON THE SEASONAL AND DAY-TO-DAY VARIATIONS 


IN THE F2 REGION 


INTRODUCTION 

Wide day-to-day variability in the F-region electron concentration have 
recently been reported, from the Thomson Scatter measurements (Mahajan 1967), 
Mahajan et al., 1968), as well as from the Doppler fading of the radio transmis- 
sions from the geostationary satellites (Garriott et al., 1965; Klobucher and 
Whitney, 1966; Titleridge, 1966; Jones, 1968). These variations can either be 
interpreted in terms of changes in the neutral particle concentration or as 
changes in the electron drift velocity originating from neutral atmospheric winds 
(King and Kohl, 1965; Kohl and King, 1967; Rishbeth, 1967) and/or from electric 
fields (Chapman and Bartels, 1940; Kato, 1956; Maeda and Kato, 1966; Fejer, 
1965). In this paper some further results on this variability are presented. 
Theoretical estimates of the changes in the neutral particle concentration and the 
drift velocity required to reproduce the observed variations are made and dis- 
cussed in terms of our current knowledge on the neutral atmosphere and the 
atmospheric winds . 

OBSERVATIONS AND DATA 

The data used are essentially the same as employed by Mahajan et al. (1968) 
for the electron content studies over Arecibo (18.4°N, 67°W) and cover the period 
Sept. 1965 through Aug. 1966. These consist of the back scatter power measure- 
ments in the range 100 - 700 Km. Pulse transmissions of 500 /as were used and 
sampled at intervals of 20 / a s with the digital autocorrelation technique developed 
by Perkins and Wand (1965). The electron concentration (N h ) for any height h is 
obtained from the measured power (P h ) as: 
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A X P h x h 2 
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where A is a constant involving primarily the parameters of the transmitting and 
receiving system, and (T e /T. ) is the electron to ion temperature ratio at the 
height h. The ratio T e /T. is deduced from the power spectra recorded simul- 
taneously in the height range 250 - 550 Km. Outside this range T e /T i is 
obtained by linear interpolation by assuming this ratio to be unity at 100 and 
1000 Km. As the constant A can not be obtained accurately, it is eliminated by 
forming a normalized electron concentration profile. The absolute scale of the 
profile is provided by scaling NmF2 from the simultaneously recorded C-3 
ionograms. 

Figure 1 is a plot of a number of electron density profiles on various days 
for about the same local time (11-13 hours). Significant change in the electron 
concentration from one observation to the next can be noted. This variability is 
examined in Figure 2 by plotting NmF2 against the observed peak height. We 
have tried to distinguish the various seasons by assigning different symbols for 
summer, winter and equinox. A strong overall positive correlation between 
hmF2 and NmF2 is evident, with some seasonal change in this correlation. The 
winter NmF2 and hmF2 values are consistently lower than the summer and 
equinox. In the sections to follow we shall try to reproduce the observed hmF2 - 
NnF2 correlation theoretically by changing the neutral concentration and then 
the wind velocity. 
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THEORY 

We shall consider the F2 region to consist of oxygen ions within a neutral 
atmosphere of predominantly oxygen atoms. The oxygen ion concentration in 
the F2 region is controlled by photo ionization with production rate Q, charge 
transfer to molecules with loss rate L[0 + ] and transport of 0 + with a velocity 
v Q+ . Assuming that the transport of ionization is only important parallel to the 
magnetic field, the steady state continuity equation has the form: 

Q - L [0 + ] - sin I |[0 + ] v Q+ J = 0 (2) 

where I is the dip angle and z the variable altitude in the plane stratified 
ionosphere. 

If one assumes that the plasma temperatures are constant with altitude, the 
motion equation has the form 

<V o [° + ] ^ ( v o+ ■ v o) = “ sin 1 ( kT 3 ° h " + m 0 g [0 + ]) (3) 

where: 

^o+o = drag coefficient of 0 + in 0 

v Q = component of the wind velocity of 0 parallel to the magnetic field 
T = T e + T. , the electron and ion temperature 
m 0 = mass of 0 and 

g = gravitational acceleration. 

Substituting (3) in (2), one gets: 


3 


( 4 ) 


sin 2 I 


B 
B z 


1 

Vo M 



+ L [0 + ] 



n sin 1 bT ( v ° [0+] ) " Q 


When v 0 = 0, Equation (4) can be solved analytically by its transformation 
to a Bessel differential equation (see e.g. Yonezawa, 1958; Geisler and Bowhill, 
1965) if one assumes that the functions Q and L decrease exponentially as: 

z - z 0 

Q = Q 0 e - "V — (5) 

z - Z 0 

L = L 0 e - — (6) 


Here q is the scale height of O, the major ionizeable constituent and i the 
scale height of N 2 or 0 2 to which 0 + transfers its charge. This implies that the 
gas temperature is considered constant with altitude and that the attenuation of 
the ionizing solar UV radiation is neglected. 


For the case v 0 f 0, Equation (4) can only be analytically solved if v 0 varies 
in a fashion such that: 

K 

V ° ~ [0] ' ^ * s a constant 


Under this condition a constant force “ s - n j , in addition to the gravitational 
force, would be applied to the ions. As this situation is not expected to exist, 
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Equation (4) has to be solved numerically. We have accomplished this by 
a perturbation method as summarized below . 

The term sin I [0 + ] v 0 ^ was considered as a perturbation in Eq 
(4), and it was applied in a successive iteration, such that in each iteration 
its effect was small. This method can be described in analytical form thrc 
system of equations of the form as Equation (4) but with variable paramete 

^ (0 5 < 1 ) 

e([0*] 0 ) = - Q 

">([0*],) = -Q + e, *inl£ (to*], v„) 
d([ 0*] 2 ) = - <5+ e 2 sin I -4- /[O*] , v 0 ) 


D ( [0,) k) = - 0 + £ k sin v 0 j 

The first equation in (7) for e 0 = 0 corresponds to the case v 0 = 0 and can be 
solved as previously described to give [0 + ] 0 . With this zero order solution, 
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the second equation in (7) is defined. This equation is then solved by expanding 
£j sin I -y^-^[0 + ] v 0 j in a series of exponential functions (two in our case) for 
each of which an analytical solution can be obtained. Depending on the magnitude 
of v Q , the parameter is such chosen that the difference to the previous solu- 
tion [O + ] 0 is small compared to [O + ] 0 . In this way the validity of the perturbation 
method is fulfilled. Through continuous application of this procedure and simul- 
taneously increasing e k as 

C 0 < € 1 < e 2 < € 3 • • ' ‘ - 1 ' 

this parameter approaches 1 as the Equation (7) approaches the Equation (4) to 
be solved. After this point, further iterations employ the correct form of 
Equation (4) and lead to a convergence of successive solutions [0 + ] k which 
thus represents a solution of Equation (4). 

RESULTS 

We shall, for simplicity, assume a constant ionizing flux. While changes in 
the flux proportionally change NmF2, these do not affect the hmF2 for the equi- 
librium profiles (see Equation (4)). Since 0 and N 2 are the major neutral 
constituents controlling the F-region ionization, we shall examine the equi- 
librium profiles by changing the number density of these particles alone. As 
0 2 has been measured to be less than N 2 by a factor of 15 at 150 Km by Hall et 
al (1967), and since the reaction rate of 0 + to 0 2 is faster thanO + to N 2 only 
by a factor of 10, we shall not consider any effects of electron loss through 0 2 . 
Taking z Q to be 200 Km, we use the following constants 

1. Photon flux = 3 x 10 10 cm -2 s“ 1 

2. Ionization cross section of 0 = 10~ 17 cm 2 
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3. 

Rate coefficient of 0 + to N 2 

= 2 x 

10 -12 s" 1 

4. 

Drag coefficient of 0 + in 0 

II 

CO 

X 

10 -32 gm cm -3 s" 1 

5. 

Sin I (Arecibo, dip = 52°) 

= 0.8 



An electron to ion temperature ratio of about 1-5 (T = T n = 900° K, T e = 1500° K) 
is assumed throughout the F-region. This sufficiently accounts for the thermal 
nonequality, because the F-region electron concentration is not very sensitive to 
T e changes (Thomas, 1966), although wide day-to-day variability in the electron 
temperature can exist (see e.g., Mahajan 1967) . 

We have changed the 0 and N 2 concentrations in systematic steps. First the 
number density of 0 is varied for a fixed ratio of [0 ] /[ N 2 ] and later [0 ] /[N 2 ] 
is varied for the previously used [0] . The results of these calculations are 

summarized in Figure 3 where a theoretical grid of NmF2 - hmF2 relationship 

[ 0 ] 

has been generated for various [0 ] and -p— p values. Figure 3 indicates that to 
reproduce the average variability observed in Figure 2, one would need more 
than a factor of 5 change in the O and N 2 concentration with a fixed [0]/[N 2 ] 
ratio of about 2-0. 

The effects of the wind on the equilibrium profiles are seen in Figure 4. 
Meridional winds with velocities of +15, 0, -15 and -100 m/sec parallel to the 
field lines have been used. The concentrations of 0 and N 2 employed are the 
experimentally measured values by Hall et al (1967) during March 1966 and are 
as follows: 

[0] = 4.4 x 10 9 cm -3 

[N 2 ] - 2.3 x 10 9 cm -3 

As can be seen that the winds have a very strong effect on the electron concen- 
tration profiles. These effects are similar to the vertical drift effects reported 
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by Yonezawa (1958) and Rishbeth and Baron (1960). One can immediately see 
from Figure 4 that changes in wind velocity can reasonably account for the 
average variability observed in the electron concentration. It is also clear that 
one would need equatorialwards winds (positive velocities) during summer and 
polewards winds (negative velocities) during winter at Arecibo to explain the 
observed seasonal changes in NmF2 and hmF2. 

Another result evident from Figure 4 is that the effect of poleward winds on 
the electron concentration profiles is relatively smaller than that of the equa- 
torialward winds. This is expected to be so since high poleward winds bring the 
F-layer peak close to the chemical equilibrium region. 

DISCUSSION 

It is clear from Figure 3 that to explain the observed winter to summer 
change in NmF2 and hmF2, one would need a factor of 5 change in the 0 and N 2 
concentrations and thus in the atmospheric density. As a matter of fact no sig- 
nificant seasonal changes in the atmospheric density from the satellite drag data 
have been reported (see Priester et al, 1967 for a review on satellite drag re- 
sults). Also the ion temperature measurements at 250 km (T. = T n at 250 km, 
Nisbet 1967) during these observations indicate neither any appreciable seasonal 
change in T. nor any systematic dependence of hmF2 (and/or NmF2) on T. . Thus 
it is difficult to invoke changes in the neutral concentration as the cause of the 
F-region electron concentration variations. 

As we have seen, the atmospheric winds can reasonably account for the ob- 
served variability. Geisler (1966) has calculated the model wind systems for the 
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F-region for mid latitude stations. He has found poleward winds during the day- 
time with velocities higher in winter than summer. As the diurnal bulge in atmos- 
pheric density (and thus in pressure) migrates in latitude in phase with the sub- 
solar point (Jacchia, 1965; Jacchia and Slowey, 1968), the bulge is expected to be 
noi'th of Arecibo during most of the summer. This would give rise to equatorial- 
ward winds — a conclusion in conformity with our results in Figure 4. Due to 
the location of the diurnal bulge to the south of Arecibo in winter and equinox, one 
would then expect poleward winds, which is again consistent with our results. 

The cause of the changes in the wind speed during the same season is, however, 
not known. It does seem evident that small changes in the wind velocity during 
summer and equinox (positive velocities or low negative velocities) would result 
in relatively greater change in the electron concentrations than during winter. 

The wide day-to-day variability (about a factor of three) in the electron content 
data of Klobuchar and Whitney (1966) during the summer of 1965 lends support to 
this prediction. 

We have implicitly neglected any effects of the electric fields in our calcula- 
tions, since we have assumed the ionization to move along field lines. Electric 
fields, if present, would have a similar effect on the ionization transport as the 
meridional winds. However, in the absence of any measurements, it is rather 
difficult to test their effects on the F2 region. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Electron concentration profiles from Thomson scatter measurements 
at Arecibo, Puerto Rico. All profiles correspond to about the same 
local time. The day-to-day and seasonal changes in the concentration 
can be noted. 

A plot of NmF2 against hmF2 for observations on various days during 
the period Sept. 1965 to Aug. 1966. Times between 11-13 hours have 
been used. 

Theoretical relationship between hmF2 and NmF2 for various values 
of [0 ] and [0 ] /[N 2 ] . The values of the other parameters used in 
the continuity equation are given in the text. 

Theoretical electron density profiles showing the wind effects. The 
0 and N 2 concentrations used have been taken from measurements 
by Hall et al (1967) and very nearly correspond to the same solar 
activity conditions as our N-h data. The velocities are the compo- 
nent along the field lines. The horizontal velocities will be twice 

the numbers ( i.e. — : — ^ 7 ) shown in the diagram. 
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Figure 1. Electron concentration profiles from Thomson scatter measurements at Arecibo, Puerto 
Rico. All profiles correspond to about the same local time. The day-to-day and seasonal changes 
in the concentration can be noted. 
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Figure 2. A plot of NmF2 against hmF2 for observations on various days during the period 
Sept. 1965 to Aug. 1966. Times between 11-13 hours have been used. 
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Figure 3. Theoret : cal relationship between hmF2 and NmF2 for various values of [0] and 
[0]/[N 2 ]. The values of the other parameters used in the continuity equation are given in the 
text. 
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ELECTRON CONCENTRATION (/ CC) 

Figure 4. Theoretical electron density profiles showing the wind effects. The 0 and N 2 concen- 
trations used have been taken from measurements by Hall et al (1967) and very nearly correspond 
to the same solar activity conditions as our N— h data. The velocities are the component along 

i 1 

the field lines. The horizontal velocities will be twice the numbers li.e. 'f T 
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